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High  resolution  simulations  of  Arctic  sea  ice,  1979-1993 

Wieslaw  Maslowski  & William  H.  Lipscomb 

"OOIARIS^' 

To  evaluate  improvements  in  modelling  Arctic  sea  ice,  we  compare 
results  from  two  regional  models  at  1/12°  horizontal  resolution.  The  first 
is  a coupled  ice-ocean  model  of  the  Arctic  Ocean,  consisting  of  an  ocean 
model  (adapted  from  the  Parallel  Ocean  Program,  Los  Alamos  National 
Laboratory  [LANL])  and  the  “old”  sea  ice  model.  The  second  model  uses 
the  same  grid  but  consists  of  an  improved  “new”  sea  ice  model  (LANL/ 

CICE)  with  a simple  ocean  mixed  layer.  Both  models  are  forced  with 
European  Centre  for  Medium-range  Weather  Forecasts  reanalysis  data 
for  1979-1993.  A comparison  of  the  two  sea  ice  models  focuses  on  the 
winter  of  1987  to  emphasize  the  internal  ice  stress  and  to  minimize  biases 
towards  a particular  Arctic  climate  regime.  The  “new”  sea  ice  model 
gives  improved  ice  deformation  and  drift  fields.  These  improvements  are 
associated  at  least  in  part  with  the  multi-category  representation  of  the 
ice  thickness  distribution  and  more  realistic  parameterization  of  the  ice 
strength.  Long,  narrow  features  in  ice  divergence  and  shear  fields  resem- 
ble those  observed  in  SAR  imagery,  except  that  their  average  width  is 
overestimated,  possibly  due  to  insufficient  horizontal  resolution.  We  also 
compare  the  mean  sea  ice  drift  and  its  decadal  variability  in  two  “old”  sea 
ice  models  at  different  horizontal  resolutions:  18-km  and  9-km.  We  find 
no  significant  change  in  ice  drift  between  the  two  models,  except  in  areas 
of  significant  ice-ocean  interactions  due  to  more  realistic  ocean  currents 
and  water  mass  properties  in  the  9-km  model. 

W.  Maslowski,  Dept,  of  Oceanography,  Naval  Postgraduate  School,  833  Dyer  Rd..  Monterey,  CA  93943,  USA, 
maslowsk@nps.navy.mil ; IV.  H.  Lipscomb,  T-3  Fluid  Dynamics  Group,  Los  Alamos  National  Laboratory,  Los 
Alamos,  NM 87545,  USA. 


With  its  sea  ice  cover,  the  Arctic  Ocean  has  been 
difficult  to  model  in  both  global  climate  simula- 
tions and  regional  studies.  One  challenge  is  to 
represent  faithfully  the  physical  processes  spec- 
ific to  sea  ice  and  oceans  in  polar  regions.  Anoth- 
er challenge,  intimately  linked  to  the  first,  is 
to  resolve  small-scale  features  such  as  narrow 
boundary  currents  (0[100km]),  the  Rossby  radius 
of  deformation  (0[10km]),  and  the  bottom  bathy- 
metry and  land  geometry  controlling  the  phys- 
ics and  inter-basin/ inter-ocean  communications. 
Climate  change  studies  suggest  that  the  Arctic  is 
highly  sensitive  to  greenhouse  warming,  large- 
ly because  of  positive  feedbacks  associated  with 
thinning  and  retreat  of  the  insulating,  highly 


reflective  sea  ice  cover.  The  ice  models  in  these 
studies  have  typically  been  run  at  fairly  coarse 
resolution  (~3°),  using  crude  parameterizations  of 
the  thermodynamic  and  dynamic  processes  that 
determine  the  ice  thickness  and  extent.  Howev- 
er, recent  advances  in  sea  ice  modelling,  together 
with  parallel  supercomputers  for  high  resolution 
modelling,  make  it  feasible  to  improve  the  fidelity 
of  Arctic  Ocean  simulations  and  thereby  increase 
the  reliability  of  climate  change  predictions. 

In  this  paper  we  present  results  from  two 
regional  models  of  the  Arctic  Ocean  at  9-km 
(or  1/12°)  resolution.  We  briefly  describe  the  ice 
and  ocean  models,  atmospheric  forcing,  and  ini- 
tial and  boundary  conditions  used  in  this  study. 
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We  then  compare  model  results,  focusing  on  sea 
ice  concentration,  velocity,  shear  and  divergence 
fields  that  can  be  qualitatively  compared  to  sim- 
ilar sea  ice  products  available  from  high  reso- 
lution SAR  imagery,  such  as  the  RADARSAT 
Geophysical  Processor  System  (RGPS)  analy- 
ses (e.g.  Kwok  2001).  To  better  understand  model 
improvements  following  from  increased  resolu- 
tion, we  also  compare  the  mean  sea  ice  drifts  and 
anomalies  from  our  earlier  18-km  and  30  level 
coupled  ice-ocean  model  (Maslowski  et  al.  2001) 
with  the  similar  coupled  model  run  at  9-km  with 
45  ocean  levels.  We  conclude  by  discussing  pos- 
sible sources  of  model  differences  and  suggesting 
further  experiments.  Due  to  length  constraints  in 
this  publication,  only  limited  presentation  of  the 
results  is  possible. 

Model  Description 

Coupled  ice-ocean  model 

The  9-km  model  uses  a parallel  version  of 
the  basic  sea  ice  model  of  Hibler  (1979),  with 
viscous-plastic  dynamics,  simple  thermodynam- 
ics and  two  thickness  categories  (thick  ice  and 
thin  ice/open  water).  The  ice  model  is  coupled  to 
a regional  adaptation  of  the  Parallel  Ocean  Pro- 
gram (POP)  of  Los  Alamos  National  Laboratory 
(LANL).  This  ocean  model  has  a free-surface  for- 
mulation (Smith  et  al.  1992)  and  45  vertical  levels. 
The  model  domain  is  defined  in  a rotated  spher- 
ical coordinate  system  to  eliminate  the  singular- 
ity at  the  North  Pole  and  was  chosen  to  include 
all  the  northern  ice-covered  oceans.  It  extends 
from  the  North  Pacific  at  ~30°N,  across  the  cen- 
tral Arctic  and  sub-Arctic  seas,  into  the  North 
Atlantic  to  ~45°N.  The  model  uses  significantly 
improved  bathymetry  data,  including  the  2.5  km 
resolution  International  Bathymetric  Chart  of  the 
Arctic  Ocean  database  (Jakobsson  et  al,  2000), 
along  with  the  recently  upgraded  monthly  mean 
temperature  and  salinity  climatology  (PHC)  from 
the  University  of  Washington  (Steele  et  al.  2001). 
The  model  includes  run-off  from  most  major 
rivers  (the  Yukon,  Mackenzie,  Dvina,  Pechora, 
Ob,  Yenisey,  Kotuoy,  Lena,  Kolyma  and  Indigir- 
ka) as  well  as  ungauged  run-off.  Other  model  for- 
mulations are  similar  to  those  used  in  the  18-km, 
30-level  coupled  ice-ocean  model  of  the  Arctic 
Ocean  of  Maslowski  et  al.  (2000). 

We  have  so  far  completed  a 63-year  integration 
of  the  coupled  9-km  ice-ocean  model.  The  model 
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was  forced  with  daily-average  atmospheric  fields 
(downward  longwave  and  shortwave  radiation,  air 
temperature  and  density,  specific  humidity,  precip- 
itation, and  wind  velocity)  from  the  European 
Centre  for  Medium-range  Weather  Forecasts 
(ECMWF)  1979-1993  reanalysis.  Surface  fluxes 
of  turbulent  heat  and  momentum  were  computed 
using  bulk  formulas.  Following  the  27-year 
spinup,  we  performed  an  additional  21-year  run 
using  the  repeated  1979  ECMWF  annual  cycle 
for  the  first  6 years  and  1979-1981  inter-annual 
fields  for  the  last  15  years.  In  this  way  we  forced 
the  sea  ice  and  ocean  toward  conditions  of  the  late 
1970s  and  early  1980s,  which  were  used  to  start 
the  1979-1993  integration. 

CICE  model 

The  “new”  sea  ice  model,  CICE  version  3.0 
(Hunke  & Lipscomb  2001),  was  developed  at 
LANL  for  climate  studies  and  differs  from  the 
old  Hibler  model  in  several  respects.  The  new 
ice  model  incorporates  the  latest  version  of  the 
elastic-viscous-plastic  (EVP)  dynamics  of  Hunke 
& Dukowicz  (1997),  which  reduces  to  the  viscous- 
plastic  model  at  long  time  scales  but  uses  elas- 
tic waves  to  improve  the  response  on  short  time 
scales  and  increase  computational  efficiency.  The 
wind  stress  depends  not  only  on  the  wind  veloc- 
ity, but  also  on  a drag  coefficient  that  varies  with 
the  boundary  layer  stability  following  Bryan  et  al. 
(1996).  (The  old  ice  model  uses  a constant  drag 
coefficient.)  Vertical  heat  transfer  and  growth  and 
melting  of  ice  and  snow  are  computed  using  the 
multi-layer  energy-conserving  thermodynamic 
model  of  Bitz  & Lipscomb  (1999),  which  includes 
a temperature-dependent  heat  capacity  to  mimic 
brine  pocket  effects.  The  ice  is  transported  hori- 
zontally using  the  incremental  remapping  scheme 
of  Lipscomb  & Hunke  (unpubl.  ms),  which  is 
second-order  accurate  in  space  except  where  the 
accuracy  is  reduced  locally  to  preserve  monoto- 
nicity. The  old  model  uses  a second-order  centred 
difference  scheme  that  requires  a diffusion  term 
to  maintain  numerical  stability. 

Perhaps  the  most  important  difference  is 
that  the  new  model  has  multiple  ice  thickness 
categories  (along  with  open  water).  Ice  is  trans- 
ferred between  categories  as  it  grows  and  melts 
using  the  thickness-remapping  scheme  of  Lips- 
comb (2001).  Several  studies  (e.g.  Maykut  1982; 
Bitz  et  al.  2001)  have  shown  that  the  ice  strength, 
heat  exchange  with  the  atmosphere,  and  other 
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properties  are  sensitive  to  the  ice  thickness 
distribution,  especially  to  the  amount  of  thin  first- 
year  ice.  With  multiple  categories,  sea  ice  ridging 
can  be  represented  as  a process  that  moves  ice 
from  thin  categories  to  thicker  categories,  as 
suggested  by  Thorndike  et  al.  (1975)  and  Hibler 
(1980).  In  the  new  model  the  ice  strength  P is 
proportional  to  the  increase  in  potential  energy 
during  ridging,  following  Rothrock  (1975).  In  the 
old  model  the  strength  is  represented  as: 

P = P*  h exp  [-C(l -/!)], 

where  h is  the  mean  ice  thickness,  A is  the  open 
water  fraction,  and  P*  and  C are  empirical 
constants  that  are  tuned  to  give  reasonable 
ice  strengths.  Following  Hibler  (1979),  we  set 
/>*=27500  N/m2  and  C=20.  In  both  parameter- 
izations,  the  ice  weakens  as  the  open  water 
fraction  increases.  In  the  old  model,  however,  the 
ice  strength  depends  on  the  mean  ice  thickness, 
whereas  in  the  new  model,  the  ice  strength 
depends  on  the  amount  and  thickness  of  thin  ice. 

We  have  completed  a 20-year  integration  of 
CICE  on  the  9-krn  grid.  The  run  consists  of  a 10- 
year  spinup  and  a 10-year  interannual  run  forced 
with  the  ECMWF  1979-1988  reanalysis  data. 
The  model  was  run  with  five  ice  thickness  cat- 
egories, each  having  four  vertical  ice  layers  and 
one  snow  layer,  and  other  standard  parameters 
as  described  in  Hunke  & Lipscomb  (2001).  The 
ice  was  coupled  not  to  a full  ocean  model,  but  to 
a simple  ocean  mixed  layer  that  can  absorb  heat 
and  grow  ice  in  leads.  The  ice  receives  heat  from 
the  ocean  at  its  bottom  surface,  although  not  as 
much  as  in  coupled  models  where  there  can  be  an 
upward  heat  flux  from  beneath  the  mixed  layer. 
Also,  the  ocean  velocity  is  set  to  zero  in  the  water 
stress  term  of  the  sea  ice  momentum  equation. 
We  emphasize  that  disparities  in  sea  ice  behav- 
iour between  the  old  and  new  models  could  arise 
not  only  from  differences  in  the  ice  models  them- 
selves, but  also  from  different  coupling  to  the 
ocean  and  atmosphere.  In  the  near  future  we 
plan  to  couple  CICE  with  POP  so  as  to  make 
more  detailed  comparisons  of  the  old  and  new  ice 
models.  The  stand-alone  integration,  however, 
already  reveals  significant  differences  between 
the  two  ice  models. 

Results 

To  differentiate  the  effects  of  a new  ice  model 


from  those  of  increased  spatial  resolution,  we  first 
compare  sea  ice  drift  in  two  coupled  ice-ocean 
models  of  the  Arctic  Ocean,  both  using  the  “old” 
ice  model  but  at  different  resolution.  The  coarser 
model,  with  an  18-km  grid  and  30  vertical  levels, 
extends  from  (closed)  Bering  Strait  to  ~45°N  in 
the  North  Atlantic.  Maslowski  et  al.  (2001)  have 
shown  that  the  mean  sea  ice  drift  in  this  model 
changes  its  circulation  over  a 10-year  period 
between  1981-82  and  1991-92.  In  the  1980s 
the  large-scale  anticyclonic  Beaufort  Gyre  pre- 
vails, with  the  Transpolar  Drift  positioned  along 
the  Lomonosov  Ridge  (Fig.  2a  in  Maslowski  et 
al.  2001).  In  the  early  1990s  the  Beaufort  Gyre 
shrinks  to  occupy  only  the  Beaufort  Sea,  and  the 
Transpolar  Drift  shifts  eastward,  aligning  itself 
along  the  Alpha  and  Mendeleyev  ridges  (Plate 
3b  in  Maslowski  et  al.  2000).  The  difference 
between  the  two  2-year  mean  ice  velocity  fields 
illustrates  a basin-wide  cyclonic  trend  (Fig.  la  in 
Maslowski  et  al.  2001).  Figures  la-c  show  equiv- 
alent 9-km  model  results  for  the  2-year  mean  sea 
ice  drift  of  1981-82  and  1991-92  and  the  differ- 
ence between  them.  The  general  large-scale  ice 
circulations  and  the  decadal  shift  between  the 
two  2-year  means  are  similar  to  results  from  the 
previous  model. 

A significant  difference  exists  in  the  Chukchi 
and  southern  Beaufort  Seas  in  1991-92,  where 
the  flow  is  anticyclonic  in  the  18-km  model  but 
cyclonic  in  the  9-km  model  (Fig.  lb).  The  18-km 
model  shows  a much  stronger  southward  flow 
over  the  Chukchi  Sea  towards  Bering  Strait  in 
the  1980s  compared  to  the  9-km  model  (Fig.  la). 
Since  both  models  use  the  same  ECMWF  forcing, 
we  attribute  these  differences  to  the  upper  ocean 
dynamics,  which  are  strongly  affected  by  the  flow 
of  Pacific  Water  northward  through  Bering  Strait, 
then  eastward  along  the  Chukchi  and  Beau- 
fort shelves.  The  9-km  model  predicts  net  mean 
northward  flow  through  the  open  Bering  Strait 
on  the  order  of  0.7 -0.8  Sv  (McClean  et  al.  2001), 
in  good  agreement  with  observations  (Roach  et 
al.  1995).  The  9-km  model  also  yields  more  real- 
istic shelf  and  boundary  ocean  currents  in  this 
region,  which  through  ice-ocean  drag  affect  the 
ice  motion.  From  this  analysis  we  conclude  that 
spatial  resolution  does  not  significantly  improve 
results  from  the  old  sea  ice  model,  except  indi- 
rectly through  coupling  of  a more  realistic  ocean 
to  sea  ice. 

The  old  model  has  a reasonable  ice  thickness 
distribution  (Fig.  Id),  but  we  find  problems  in 
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Fig.  I.  The  mean  sea  ice  thickness  distribution  (m)  and  drift  (m/s)  for  (a)  1981-82,  (b)  1991-92,  (c)  the  difference  between 
1991-92  and  1981-82,  and  (d)  the  1979-1993  mean.  Thick  contour  lines  show  the  ice  thickness  every  1 m between  1 and  5 m. 
Thin  contours  in  (d)  are  every  1 m between  0.5  and  4.5  m 


the  velocity  and  deformation  fields.  Figure  2 
shows  snapshots  of  ice  concentration  and  drift 
(2a,  b),  divergence  (2c,  d),  and  shear  (2e,  f)  from 
the  old  ice  model.  These  snapshots  are  taken  six 
days  apart  to  show  the  response  to  a synoptic- 
scale  change  in  atmospheric  forcing.  We  chose 
results  from  mid-winter  1987  to  maximize  the 
influence  of  the  internal  ice  stress  and  to  avoid 
a bias  toward  either  the  cyclonic  or  anticyclonic 
atmospheric  regime.  The  left  and  right  columns 
show  fields  for  20  and  26  February  1987,  respec- 
tively. Compared  to  buoy-measured  ice  speeds 
(http://iabp.  apl.washington.edu/),  which  are  typ- 
ically in  the  range  of  0-10  cm/s  in  the  central 
Arctic,  the  ice  speeds  are  too  small.  Large  por- 
tions of  the  winter  ice  pack  are  almost  complete- 
ly still,  even  with  ice-ocean  drag  present.  The 
shear  and  divergence  fields  lack  the  long,  narrow 
features  seen  in  RGPS  data.  Instead,  most  of  the 


action  takes  place  in  large  diffuse  patches  of  high 
positive  divergence  and  shear. 

Figure  3 shows  corresponding  fields  for  the 
new  ice  model.  In  Fig.  3a  and  b,  ice  speeds  in 
most  of  the  Arctic  basin  are  an  order  of  magni- 
tude larger  than  in  the  old  model,  in  good  agree- 
ment with  buoy  measurements.  The  anticyclonic 
Beaufort  Gyre  is  present  in  buoy  data  for  Febru- 
ary 1987,  though  not  as  strong  as  in  the  plot  from 
20  February.  (The  model  skill  validated  against 
buoy  motions  on  a particular  day  depends  strong- 
ly on  the  quality  of  the  smoothed,  coarse  resolu- 
tion atmospheric  forcing.)  A significant  change  in 
large-scale  ice  drift  over  6 days  suggests  a strong 
dependence  on  synoptic-scale  winds.  The  model 
ice  concentration  in  the  central  Arctic  is  close 
to  100%  as  expected.  The  divergence  and  shear 
fields  show  long,  narrow  features  with  sharp  grad- 
ients, which  are  qualitatively  similar  to  RGPS 
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Fig.  2.  Snapshots  of  (a,  b)  the 
total  sea  ice  concentration 
(%)  and  drift  (m/s),  (c,  d)  sea 
ice  divergence  (%/day),  and 
(e,  f)  shear  (%/day)  on  (left) 

20  February  and  (right)  26 
February  1987  from  the  coupled 
ice-ocean  model  with  the  “old” 
sea  ice  model.  Unit  drift  vector 
is  0.1  m/s. 


observations  (Kwok  2001).  Upon  closer  com- 
parison of  high  resolution  SAR  imagery  to  model 
deformation  features,  we  find  the  model  features 
to  be  wider  and  less  frequent  in  space  and  time. 
We  hypothesize  that  the  model  spatial  resolution 
is  insufficient  to  resolve  the  finer  features  meas- 
ured with  RGPS.  Nevertheless,  the  “new”  sea  ice 
model  shows  dramatic  improvements  compared 
to  the  old  model.  We  offer  possible  explanations 
in  the  next  section. 

Discussion  and  conclusions 

The  18-km  and  9-km  “old”  sea  ice  models  have 
similar  skill  in  representing  the  climate  regime 
shift  in  the  Arctic  Ocean  from  the  early  1980s 
to  the  early  1990s.  The  similarity  in  ice  drift 
between  the  two  models  suggests  that  increas- 


ing the  model  resolution  has  little  effect  on  the  ice 
strength.  (One  might  have  predicted  that  the  ice 
strength  and  resulting  drift  would  change  due  to 
increased  resolution  of  fine-scale  features  such  as 
leads.)  There  are  some  improvements  in  the  repre- 
sentation of  ice  edge  position,  polynyas,  concent- 
ration and  thickness  compared  to  lower  resolution 
sea  ice  models  (not  shown).  However,  many  of  the 
improvements,  such  as  the  ice  edge  position  in  the 
Nordic  and  Bering/Chukchi  seas,  are  influenced 
by  more  realistic  modelling  of  the  ocean  currents 
and  hydrography  at  increased  resolution. 

The  old  ice  model  gives  ice  speeds  that  are 
too  small,  despite  having  produced  reasonable 
velocity  fields  in  other  studies,  including  Hibler 
(1979).  In  this  study  the  mean  ice  thickness 
is  3-4  m in  the  central  Arctic,  giving  an  ice 
strength  P of  the  order  of  105  N/m  during  winter 
when  the  open  water. fraction  is  small.  In  Hibler 
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Fig.  3.  Snapshots  of  (a,  b)  the 
total  sea  ice  concentration  (%) 
and  drift  (m/s),  (c,  d)  sea  ice 
divergence  (%/ day),  and  (e,  0 
shear  (%/day)  on  (left)  20  Febru- 
ary and  (right)  26  February  1987 
from  the  stand-alone  “new”  sea 
ice  model.  Unit  drift  vector  is 
0.1  m/s. 


(1979)  most  of  the  ice  was  thinner  than  3 m, 
giving  much  lower  ice  strengths.  As  shown  by 
Steele  et  al.  (1997),  the  ice  motion  field  is  a highly 
nonlinear  function  of  ice  strength.  These  authors 
found  that  ice  moves  at  reasonable  speeds  when 
the  average  strength  P is  about  3 * 104  N/m,  but 
rapidly  locks  up  as  the  strength  is  increased.  The 
mean  ice  thickness  in  our  study  is  somewhat 
larger  than  observed,  possibly  due  to  biases  in  the 
atmospheric  forcing. 

To  test  the  sensitivity  of  the  velocity  and  defor- 
mation fields  to  the  mean  ice  thickness,  we 
repeated  the  simulation  for  the  first  two  months 
of  1987 — but  with  the  ice  thickness  cut  in  half  on 
1 January,  giving  a thickness  of  less  than  2 m in 
most  of  the  Arctic.  The  resulting  fields  for  26  Feb- 
ruary are  shown  in  Fig.  4a,  c and  e,  which  should 
be  compared  to  Fig.  2b,  d and  f.  Ice  speeds  are 
generally  much  larger,  though  still  low  on  the 
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Pacific  side  of  the  Arctic.  We  obtained  similar 
results  for  20  February  (not  shown).  We  con- 
clude that  the  ice  strength  and  resulting  drift  in 
the  old  model  are  very  sensitive  to  the  mean  ice 
thickness.  The  value  of  the  strength  parameter 
P*  could  be  lowered  to  give  more  realistic  veloc- 
ities in  the  presence  of  thicker  ice,  but  we  would 
prefer  to  avoid  such  arbitrary  tuning.  The  defor- 
mation fields  do  not  improve  significantly  when 
the  ice  is  weaker  and  flowing  faster.  These  fields 
have  more  structure  over  a larger  region  but  are 
still  very  diffuse. 

The  new  sea  ice  model  yields  more  realis- 
tic velocity  and  deformation  fields  than  the  old 
model.  Moreover,  these  fields  are  not  very  sen- 
sitive to  the  mean  ice  thickness.  We  repeated  the 
above  sensitivity  test  with  the  new  model,  halving 
the  ice  thickness  in  each  category  on  1 January 
1987,  then  transferring  ice  to  thinner  categories 
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Fig.  4.  Snapshots  of  (a,  b)  the 
total  sea  ice  concentration  (%) 
and  drift  (m/s),  (c,  d)  sea  ice 
divergence  (%/day),  and  (e,f) 
shear  (%/day)  from  sensitivity 
runs  with  reduced  sea  ice 
thickness  by  50%  from  (left) 
the  coupled  ice-ocean  model 
with  the  “old”  sea  ice  model 
and  from  (right)  the  stand  alone 
“new”  sea  ice  model  on  26 
February  1987.  Unit  drift  vector 
is  0.1  m/s. 


as  necessary.  The  resulting  fields  on  26  Febru- 
ary are  shown  in  Fig.  4b,  d and  f,  which  should 
be  compared  with  Fig.  3b,  d and  f.  Although  there 
are  small  local  differences,  the  overall  similari- 
ty of  the  two  sets  of  fields  is  striking.  This  is  also 
true  for  20  February  (not  shown).  The  similarity 
can  be  explained  as  follows.  In  the  new  ice  model 
the  ice  strength  depends  not  on  the  mean  ice 
thickness  h but  on  the  thickness  of  ridging  ice  h,. 
(More  precisely,  the  ice  strength  is  proportional 
to  the  increase  in  potential  energy  during  ridging, 
which  itself  is  roughly  proportional  to  the  square 
of  the  mean  ridge  thickness.  The  mean  ridge 
thickness  in  turn  is  proportional  to  the  square 
root  of  h,  in  each  ridging  category,  following 
Hibler  [1980].)  Only  ice  in  the  thinnest  15%  of 
the  thickness  distribution — mainly  new  first-year 
ice — contributes  to  ridging.  Although  the  ice 
thickness  is  halved  initially,  within  a few  weeks 


the  thin  ice  categories  have  nearly  the  same  ice 
area  fraction  and  thickness  as  in  the  control  run. 
As  a result,  the  ice  strength  is  almost  unchanged. 

If  it  is  true  that  the  strength  of  sea  ice  depends 
primarily  on  the  amount  and  thickness  of  thin  ice, 
then  the  behaviour  of  the  new,  multi-category  ice 
model  is  much  more  realistic  than  that  of  the  old, 
two-category  ice  model.  Also,  the  new  model  is 
less  likely  to  require  special  tuning  for  different 
atmospheric  and  ocean  forcing.  Imagine,  for 
example,  a climate  change  that  increases  summer 
melting  without  altering  winter  growth.  The 
resulting  ice  would  be  thinner  on  average,  but  its 
winter  ice  strength  would  be  about  the  same.  The 
new  model,  unlike  the  old  model,  could  capture 
this  change  accurately. 

One  suspected  bias  in  the  ice  thickness  fields  of 
the  new  model  is  the  build-up  of  thick  ice  along 
the  Siberian  coast  (not  shown),  instead  of  along 
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the  Canadian  archipelago  and  Greenland  as  in  the 
old  model  (Fig.  1)  and  in  submarine  sonar  obser- 
vations (e.g.  Bourke  & Garrett  1987).  (We  note, 
however,  that  relatively  few  ice  thickness  obser- 
vations are  available  along  the  Siberian  coast.) 
This  bias  may  result,  in  part,  from  setting  the 
ocean  velocity  to  zero  in  the  momentum  equa- 
tion. The  winds  were  primarily  toward  Siberia 
in  the  months  before  February  1987,  and  ice  in 
the  uncoupled  new  model  probably  responds  too 
strongly  to  winds.  The  spatial  pattern  of  ice  thick- 
ness will  likely  improve  when  the  ice  model  is 
coupled  to  an  ocean  model.  We  will  then  be  able 
to  make  a more  definitive  comparison  between 
the  ice  concentration,  thickness  and  velocity 
fields  of  the  old  and  new  models. 

The  most  encouraging  result  from  the  new 
model  is  the  high  degree  of  realism  in  the  ice 
deformation  fields.  The  combination  of  high  res- 
olution with  EVP  dynamics,  multiple  ice  thick- 
ness categories,  and  realistic  ice  strength  pro- 
duces shear  and  divergence  fields  similar  to 
those  derived  from  RGPS  data.  Despite  known 
simplifications  in  the  EVP  model — for  instance, 
the  assumptions  of  isotropy  and  a normal  flow 
rule — the  rheology  succeeds  in  reproducing  fea- 
tures that  extend  hundreds  of  kilometres  across 
the  Arctic  basin  and  vary  on  synoptic  time  scales. 
Although  multi-decade  integrations  at  higher 
resolution  are  not  computationally  feasible  at 
present,  we  plan  to  repeat  these  experiments  for 
shorter  time  periods  or  over  a smaller  region,  to 
see  whether  deformation  features  become  more 
realistic  as  the  grid  scale  approaches  the  size  of 
individual  floes. 
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